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Abstract: We study domain-walls and bubble nucleation in a non-relativistic vector 
field theory with different longitudinal and transverse speeds of sound. We describe 
analytical and numerical methods to calculate the orientation dependent domain-wall 
tension, <j(8). We then use this tension to calculate the critical bubble shape. The 
longitudinally oriented domain- wall tends to be the heaviest, and sometime suffers an 
instability. It can spontaneously break into zigzag segments. In this case, the critical 
bubble develops kinks, and its energy, and therefore the tunneling rate, scales with the 
sound speeds very differently than what would be expected for a smooth bubble. 



* ali@phys.columbia.edu 
t xx2146@columbia.edu 

* isheng.yang@gmail.com 



Contents 



1. Introduction and Outline 



1 



2. Orientation Dependence 

2.1 Two Dimensions 



3 

4 
6 



2.2 Flat Wall Instability 



3. Bubble Shape 

3.1 Smooth Bubbles 

3.2 Kinky Bubbles 



8 
11 



12 



4. Conclusion 



14 



A. Examples of different potentials. 



17 



1. Introduction and Outline 

The study of first-order phase transitions is a fascinating subject that appears in many 
branches of physics. The standard picture is to nucleate (thermally or quantum me- 
chanically) a bubble in a homogeneous background of the false vacuum. The bubble 
interior is in the true vacuum, and it is surrounded by domain-walls — the minimal 
energy field interpolation between the false and true vacua. 

In this paper we will focus on thermal nucleation, where the critical bubble is the 
lowest saddle point of the energy barrier. For thermal tunneling, the time variable is 
not that important and focusing on theories which have vectors transforming under the 
spatial rotation can shed light on many of the subtleties. In the simplest example, a 
scalar field theory, one can show that the critical bubble must have SO(N) symmetry in 
N dimensional space [1,2]. This leads to the commonly used estimate for the tunneling 
rate T, 
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Here the critical bubble energy E s is determined by the domain-wall tension a and the 
energy difference AV, assuming a spherical bubble. The tension is given by the path 
in the field space which minimizes that integral. The scaling with the Fermi velocity 
Vf follows from the equation of motion. 

In this paper we will generalize the theory to include vector fields. Our motiva- 
tion comes form condensed matter systems like liquid crystals, Helium 3 and Langmuir 
monolayers [4-11]. The simplest vector fields to imagine are the non-relativistic vector 
fields transforming under the spatial rotation group. In (n + l)-dimensional space- 
time, these vectors have n components. We study the transitions between two discrete 
minima <fi± of a field with the Lagrangian 



In order to make the energy bounded from below, we need cl > ct ■ When cl 7^ cr, 
the potential can minimally break the spatial SO(N) symmetry. We will focus on the 
case with minimal breaking. For any field configuration that involves these two vacua, 

— * — * 

at least (</>+ — </>_) is a special direction that specifies the longitudinal wall and breaks 
the symmetry down to SO(N — 1). 

In Sec. 2, we study planar domain-walls. Due to the broken symmetry, the domain- 
wall tension acquires an orientation dependence. We set up the general analytical 
and numerical process to determine a (9), where 9 is defined as the angle between the 
normal vector of the wall and (0 + — </>_). We demonstrate a rich behavior of a(9) 
through examples in Appendix. A. We further show that in the orientations which the 
domain-wall is heavy, it may develop an instability and spontaneously break into zigzag 
segments of lighter walls. 

In Sec. 3, we solve for the shapes of critical bubbles from cr(9). The solution has 
a simple form when the above stated instability does not occur. When it does, the 
function describing the bubble shape becomes multi-valued. We show that it still has 
a simple interpretation and describes bubbles with kinks. We then calculate how the 
deformed critical bubble modifies the transition rate. 

The technique we use to solve for the shapes is identical to that for equilibrium 
bubbles, known as the Wulff construction [3]. It has been applied to "soft matter" 
systems like liquid crystals and Langmuir monolayers [4,8-11]. Our result agrees with 
the major conclusions in the these earlier works. In Sec. 4 we will summarize a few 
concepts sharpened by our analysis, and also provide an intuitive understanding of 
when and how the tunneling rates are modified. 




(1.3) 
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2. Orientation Dependence 



The Lagrangian in Eq. (1.3) leads to the following equation of motion, 



<j>i - c^djcpi - (c| - c^didj^j 



dV 

d<pi 



(2.1) 



where it is more apparent that ct and cl correspond to the transverse and the longi- 
tudinal sound speeds. 

We want to have two isolated vacua in V . This is quite easy to achieve using the 
following potential. 



Afterward, a small a can break the degeneracy to allow a first-order phase transition. 
This is just an example to show how achievable our setup is. Our further analysis will 
either be independent of the form of the potential, or focus on examples similar but 
even simpler than Eq. (2.2). 

To study first-order phase transitions, a useful starting point is the thin-wall bubble. 
First we pretend that the two vacua are degenerate and find an interpolation between 
them, which is a domain-wall. The property of the domain-wall will then later be used 
to form a bubble of the nucleation event. 

The important property already shows up when we consider the domain-wall. Since 
the interpolation between the two vacua is a vector in the field space, it breaks the 
spatial rotational symmetry, as shown in Fig. 1. How the vector <p + continuously 
changes into can be a complicated process and clearly depends on the orientation. In 
the thin wall approximation, we can summarize the effect as an orientation-dependent 
tension <j(8). When the tension is a constant, a first-order phase transition involves 
the nucleation of a spherically symmetric bubble. So naturally in a vector field system, 
orientation dependence of a (8) can lead to a nontrivial bubble shape. Here we will 
provide the general formalism to find cr(9), and then in Sec. 3 we will use it to find the 
bubble shape. 



V(f) = T l<3 2 + ^l 4 + *{H ■ $) + b(H • 0) 2 . 



(2.2) 



The last two terms are the two lowest orders of the effect from an external field H. 
We start by considering a = 0, then b < picks a preferred direction along H. When 
b\H\ 2 + m 2 /2 < 0, we get two degenerate vacua at 




(2.3) 
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Figure 1: The blue (longer) and red (shorter) arrows represent the vector field value of two 
vacua. The thick green line is the domain- wall. From left to right, we show a longitudinal 
wall, a transverse wall, and a wall with orientation 9. The orientation is defined such that for 
a longitudinal wall 8 = 0, and for a transverse wall 9 = tt/2. 



2.1 Two Dimensions 

We will demonstrate our technique in the simplest example — a vector field in 2D. For 
a potential with two degenerate vacua <fi±, a domain- wall is a static solution to the 
equation of motion, 

dV 

-c 2 T (d 2 x + d 2 y )<p x - (cf - c 2 T )d x (d x <p x + d y (f) y ) = -— , 

dV 

-c 2 T {dl + d 2 ^ - (c 2 L - 4)d y (d x( p x + d y <p y ) = -— . (2.4) 

The boundary condition is specified by two orthogonal vectors u ■ v = 0, such that 

(tf.V)0 = O, (2.5) 
lim f(Xu) = $± . (2.6) 

A— >±oo 

Namely, the field value interpolates between the two vacua purely along the normal 
vector of the domain-wall, u. 

For simplicity, we can actually always choose u = y and instead apply a rotation 
on the potential, 

Ve((f>x,<t)y) = V(<p x cos 9 + (f) y sin 9, (f) y cos 9 - (f) x sin 9) . (2.7) 
This simplifies the equation of motion to 



-'x 

dV e 
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The solution we get here is the y oriented domain- wall in potential Vg, which is equiv- 
alent to the u oriented domain-wall in the original potential V with u ■ y = cos 6. 
The tension of the domain-wall is given by the total energy per unit x. 

1 



a{9) = / dy 



(2.9) 



Here we set V = in the vacua. It is well-known that the practical way to find the 
domain- wall solution is to numerically minimize this tension [12-14], which is what we 
do in Appendix A. 

In principle, the orientation dependence of a can be arbitrarily complicated through 
Vg. Here we would like to start from a simple, yet in some sense typical case. Imagine 
the situation where at 9 = 0, the interpolation is purely longitudinal, <p x = const. 1 
Since a rotation of tt/2 just exchanges <fi x and <p y , the interpolation will become purely 
transverse with <fi y = const. It is then easy to work out from Eq. (2.9) that 

a(0) =c L f VW\d$\ , 

J path 

<t(J) = or f VW\d$\ , (2.10) 

J path 



2 

where the two integration paths are the same, so 

<r(0) c L 



a(f) c T 



(2.11) 



Potentials given by Eq. (2.2) when m 2 > satisfy the above assumptions, so do 
the simpler potentials we use in Appendix A. They not only show a good agreement 
with Eq. (2.11), but also demonstrate an excellent fit to a naive interpolation, 

a(9) = CT(0)cos 2 # + a(^)sin 2 # , (2.12) 

in the regular range of parameters. From the symmetry of the problem, it seems natural 
to expand o \ff) as a polynomial of sin 2 9 and keep the lowest order terms. 

We also analyze two extreme choices of parameters in Appendix A. One of them 
corresponds to the following tension. 



a{9) = ^Ax(O) 2 cos 2 9 + a(|) 2 sin 2 9 . (2.13) 

It turns out that Eq. (2.12) and (2.13) are quite representative for our further anal- 
ysis. Despite their simple forms which will simplify the calculation, they can actually 
have dramatically different behaviors. 

1 Note that we talk about a particular solution, instead of imposing some symmetry on V. This 
is necessary. One might try a rotational (reflection in the 2D case) symmetry on V along the vector 
— 4>~)- That turns out to be not necessary nor sufficient to guarantee that <p x is constant. 
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2.2 Flat Wall Instability 

Now we can think about a very practical question. Since cl > or, cr(0) is most likely 
the maximum tension. Even if the boundary condition is set up to preserve the x 
translational symmetry, the minimum energy interpolation can spontaneously break 
that symmetry. In plain words, we might be able to replace the flat wall with a large 
tension cr(0) by non-flat walls with smaller tensions, and hence reduce the total energy. 

The full treatment of this problem is to remove condition (2.5) and see if a symmetry 
breaking configuration can further minimize the total energy. That is a quite involved 
numerical work which we will not pursue in this paper. We will simply demonstrate 
this possibility in the thin-wall approximation. 

The total energy is a functional of the domain- wall shape, y(x). 



Now, given a symmetric boundary condition y(xi) = y{x-i) that naively asks for a flat 
wall, we can ask two questions: 

• Is the solution with y' = — tan 9 = a stable minimum of the total energy? 

• Is there a different solution y(x) that gives the global minimum? 

In other words, is a flat wall perturbatively and non-perturbatively stable? 
Given that er(0) is a local maximum, expanding E near 6 = gives 



Thus, the perturbative stability condition is 

±_d^_ > _ 1 (r 
a(0) M* e=Q ■ { ' 

Next, the non-perturbative instability is about whether there is a 6 ^ such that 



When either or both instability exists, there will be a critical angle 6 C such that 
(j{Q c )l cos9 c is the global minimum, and the wall prefers to settle into the zig-zag 
configuration that every segment is oriented at 9 C , as shown in Fig.2. One can perform 
the same stability analysis for the other initial angles. In this paper we will focus on 




(2.14) 




(2.15) 



a (6) < cr(0) cos 9 . 



(2.17) 
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Figure 2: When the flat, longitudinal wall on the left figure is too massive, it will sponta- 
neously breaks into zigzag segments in the right figure. Although the total wall area increases, 
the reduced tension still reduces total energy. 



simple cases where the global minimum 6 C is the only local minimum. It is straight- 
forward to see that domain- walls more massive than a(8 c ) always break into zigzags, 
while the lighter walls are unaffected. 

One may question the validity of our thin- wall analysis since usually the appearance 
of a kink means thick-wall effects are involved — it cannot be infinitely sharp and one 
needs to resolve the wall to understand it. We should remind our readers that such 
concern is not important at this point. Indeed a kink resolved by thick wall analysis 
will contribute a finite term to the total energy 2 , so it is not naively energetic favorable 
to produce them. However, in our planar wall setup, the x and y direction can be 
infinitely extended. We only need a finite number of kinks to gain an arbitrarily large 
amount of energy by turning flat walls into 8 C zigzags, so our analysis is sufficient to 
determine whether it can happen. 3 When later talking about bubbles, we will see that 
only two kinks are necessary. We can tune AV such that the bubble is arbitrarily large 
and approaches the flat wall situation, so the same logic applies. 

Now take a look at the two examples for <j(6), Eq. (2.12) and (2.13). Their stability 
features are dramatically different. If the orientation dependence is given by Eq. (2.13), 
then the flat wall is always stable. On the other hand, the tension given by Eq. (2.12) 
develops a perturbative instability once 



2 This has been studied in various examples under the name "boojum" . 

3 This is just an example of the following general concept. By definition, thin wall approximation 
works if the wall is thin — relatively to other length scales in the problem. For a bubble it is compared to 
the bubble size. In the flat wall setup, every other length scale is infinite, so the thin wall approximation 
is by definition good. 



a(0) 



> 2 . 





A flat wall with 9 = would break and settle into zigzag segments with 

Or. 



sm 



/<r(o) - 2(7(f ; 



(7(0) 



sm 
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c L - 2c 2 



(2.19) 



3. Bubble Shape 



Now we break the degeneracy between the two vacua by a small amount AV such that 
the thin wall approximation is still valid. The phase transition mediated by a thermally 
nucleated bubble has the rate 



cxp 



k h T 



(3.1) 



where E s is the saddle point energy of the bubble. We can find this saddle point by 
treating E cLS cL functional of the bubble shape y(x), 



E[y(x)] = ^(surface circs. ) — AV^(volume) 



dx 



a (- tan^ 1 y') ^1 + y' 2 - y AV 



(3.2) 



where the symmetry allows us to cut the bubble into four quadrants. We will focus 
on the first quadrant where y' = — tan 9. In order to perform functional variation, the 
standard boundary condition is y(x max ) = at an undetermined x max , and y'(0) = 0. 
From Sec. 2. 2 we learned to replace y'(0) = by y'(0) = — tan# c instead. Despite that 
it is not smooth, it does eliminate the boundary variation and is indeed what we get 
from the Euler-Lagrange equation. A more formal argument is to write down E[x(y)] 
instead, for which the standard choice, x = and x' = does not exclude the kink. 
The resulting Euler-Lagrange equation is essentially the same and we can just use it. 
We will keep it simple and solve the Euler-Lagrange equation for y(x). 



AV 



d 

dx 



da 



y 



12 



a 



•A 



V 



12 



d 

dx 



The general solution can be parametrized by 9 

1 



x{6) = const. + 



a sin 9 



da 



a sin 9 



cos 9 



da 
d9 



cos 9 



(3.3) 



(3.4) 



A^ V d9 

Note that the quantity in the parenthesis is zero at both 9 = 9 C and 9 = 0. Therefore, 
solutions starting at either value will eliminate boundary variations as promised, and 
also set that integration constant to zero. We can then integrate to find 

If da 



y{0) 



AV 



a cos 9 



dO 



sin 9 



(3.5) 
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We can try to generalize this to N dimensions, such as taking x^ as the longitudinal 
direction along (</>+ — </>_) which we will denote by xl, and x\ to xn-i as the transverse 
directions. Be aware that in general the structure of the potential can further break 
the SO(N — 1) symmetry, since the exact interpolation between <fr± may still involve 
nontrivial profile of the transverse fields. That being said, we will focus on the simple 
cases with SO(N — 1) symmetry, in which a is again only a function of 9. We can 
simply write down the energy 

E[x L {x T )\ = 2S N ^2 / dx T x*~ 2 ( o^tan^ 1 x' L )yJl + x'l - AVx L ) , (3.6) 



where Sn-2 is the area for an (N — 2) unit sphere, xt = y J^i 1 xf. 
This leads to the general solution, 

Xt(9) = (a sin 9 + ^ cos ej , (3.7) 

x L (9) = (a cose- ^sinfl) . (3.8) 

Although this is a naive generalization of [3], we should take a closer look. Note 
that by symmetry, we have ^| = at and ir/2. Also for simplicity we can treat a{9) 
as a monotonically decreasing function. So we can see that x^ is positive definite, but 
there is a risk of x being negative. Since x(n/2) is still always positive definite, and 
x goes to zero exactly at 8 C , if we naively plot Eq. (3.8) we may get something like a 
wrapped candy, as in Fig.3. 

We will provide a simple argument to prove the following statement. 

Eq. (3.8) always gives the correct critical bubble profile. When the flat 
longitudinal domain- wall is stable, it works with < < iz/2. When the 
fiat longitudinal domain-wall is unstable, we should take the largest C such 
that x T (0 c ) = and use the portion C < < ir/2. In other words, cut the extra 
wrappings and keep the candy. 

First of all, the saddle point we are looking for has only one negative mode, which 
corresponds to the expansion/contraction of the bubble. This means that fluctuations of 
the wall shape should still correspond to positive modes. So locally every wall segment 
still settles to the minimum energy configuration. As shown in Sec. 2. 2, domain-walls 
with 9 > 9 C can stay, but those with 9 < 9 C cannot exist on the critical bubble profile. 
Therefore, we have to cut off the tails, and the profile necessarily includes a kink. 

Next, can the kink occur at some 9 > # c ? Picture this in 2D for better intuition, 
that is like using a smaller portion of the two shells. Also, can we take the two shells 
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Figure 3: The bubble profile given by Eq. (3.8) and (2.12) when cl > 2ct- The correct 
profile of the critical bubble is simply the middle portion. 

further apart and interpolate between them with zigzag walls? The former possibility 
is making the bubble smaller, while the later is making it bigger. Through a pictorial 
argument, we can show that they both make the total energy smaller, establishing that 
Eq. (3.19) is really the saddle point with this unique negative mode. The foundation 
of our argument is 

AVx L (9 c ) = (N-l)^ 6 j- (3 .9) 
cos U c 

which we can get from Eq. (3.8). The physical meaning is that the energy difference due 
to volume for a cylinder — an (N — 1) sphere times height xl{@ c ), is equal to the energy 
in the domain- wall that covers the (N — 1) sphere by a zigzag profile with orientation 

e c . 

Then, as shown in the left portion of Fig.4, the energy lost due to the green (shaded) 
region is equal to the contribution from the dotted domain-wall. After they cancel 
each other, the two extra triangular regions still contribute — AV^, so the total energy 
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Figure 4: The left figure visualizes Eq. (3.10), where we attempt to make a larger bubble 
by inserting true vacuum regions and extra interpolation walls. The right figure visualizes 
Eq. (3.11), where we try to make a smaller bubble by removing some part of the walls and the 
true vacuum region. Both result in smaller total energy, which shows that the kinky shape is 
indeed a saddle point. 

is indeed less. 

(zigzag wall) = A^(rectangle) < Al/(extra false vacuum region) . (3.10) 

In the right portion of Fig.4 we try to make a smaller bubble by removing the 
true vacuum region and domain-walls covered by the green (shaded) rectangle. Then 
patching the remaining two shells together as a smaller bubble, with a kink angle larger 
than 9 C . Since cr(9 c )/ cos# c is a minimum of u{9)/ cos 9, losing those wall segments 
over-compensates the energy gain even if we remove — AV of the entire green (shaded) 
rectangle, and there are even those 4 corners that we are not really removing. So the 
energy of the resulting smaller bubble is also less. 

(removed wall) > (zigzag wall) = (rectangle) > Al^(removed false vacuum region) . 

(3.11) 

3.1 Smooth Bubbles 

We can get more intuition by solving the exact bubble shape from a specific a (9). In 
the first example we will use Eq. (2.13), where no spontaneous symmetry breaking 
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should occur. Thus, we are expecting a smooth bubble. Plugging into Eq. (3.8), we get 

N-l a(n/2) 2 sm9 



x T (9) 
x l (0) 



v /a(0) 2 cos 2 ^ + ( 7(7r/2) 2 sin 2 ^ 
N-l a(0) 2 cos9 



Ay v /a(0) 2 cos 2 ^ + ( r(7r/2) 2 sin 2 e ' 
Obviously, the bubble takes the shape of an ellipsoid, 



(3.12) 
(3.13) 



S^N-1 T 2 2 



(N - l) 2 r 2 , 



where 



r 



<r(0) o-(tt/2) 



2y|d0| 



Ci Ay op av Ay Jpath 

comes from Eq. (2.10). 

It is then straightforward to calculate the saddle point energy, 



(3.14) 



(3.15) 



E„ = S. 



(N-l) 



N-l 



A? 



iV 



Ay^- 



(3.16) 



Compare this answer to the usual form people use assuming a spherical bubble, Eq. (1.2), 
the difference can be characterized by an effective Fermi velocity, 



which is a weighted geometric average of sound speeds. 
3.2 Kinky Bubbles 

Now we turn our attention to Eq. (2.12). Plugging it into Eq. (3.8), we get 
(N-l) 



(3.17) 



x T (0) 
x L (0) 



AV 
(N-l) 
AV 



2a(n/2) - a(0) 
2(j(0) - ct(tx/2) 



sin 9 - 
cos 9 



a(0) - a(n/2) 
a(0) - o(tx/2) 



sin 3 #J , (3.18) 
cos 3 #) . (3.19) 



As expected, when <r(0) < 2cr(7r/2), we still have a smooth bubble profile. When 
<j(0) > 2<t(7t/2), as proved in Sec. 3 we just use the portion 9 C < 9 < tt/2. 
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The expression of E s is quite complicated in arbitrary dimensions, so we only 
present the "realistic" dimensions. For iV = 2, we have 



E^ 2 = ^ Path 4 ^' #1; (lOc L c T -cl- elf- , for c L < 2c T , (3.20) 

= 4AV V ( 10c i c r - c L - c T ) cos v 

+ (ex + 13c t )a/(c l - 2c T )c T ^ , for c L > 2c T 



cl - c T 



(3.21) 



This is quite complicated. We should again compare it to the spherical bubble and 
think in terms of the effective Fermi velocity, especially in the limit cl 3> ct- 

fl0c L c T -c 2 L -c 2 T \ 1/2 

VF -> ( g " J > fOT C L < 2CT , 

4 

w F -> -=(c L c T ) 1/4 , for c L > c T . (3.22) 
v3 

For A^ = 3, we have 

f n=z _ . (-U v^l^l) 3 4(c| - 10c|c T + 52c L c r - 8c T ) 

(/pathV^I^Q 32c T (7c L -6c T ) v ^ 
= 47F 10 5v ^^ ' f ° r ° L>2CT - (3 ' 23) 

Similarly we have 

v F (^ c i /2 4 /2 ) 7 , for c L > c T . (3.24) 

Comparing these to Eq. (3.17), we found that the effective Fermi velocity is still a 
weighted geometric mean, but the weight on cl is always reduced by half. This is quite 
understandable since the critical bubble approaches a thin slit. A major portion of its 
domain-wall is aligned in the transverse direction. It is straightforward to show that 
this limit generalizes to A" dimensions as 

/ 1/2 AT-l/2xl/JV (n nt -\ 

v F ~ [Cl c t ) ■ ( 3 - 25 ) 
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4. Conclusion 



We studied the orientation dependence of the domain-wall tension in a vector field 
theory. We then constructed critical bubbles for thermal nucleation. The shape of 
the bubble directly depends on the tension through a simple formula, Eq. (3.8). The 
longitudinally oriented domain-wall is usually the most massive, thus it may sponta- 
neously break into zigzag segments of a critical orientation 9 C . When that happens, the 
critical bubble develops two kinks of angle 9 C , and the overall shape is still described 
by Eq. (3.8) with a careful interpretation. 4 

Our analytic and numerical study shows that the freedom to take different paths in 
a multi-dimensional field space is essential for the instability. If we choose parameters 
that reduce the number of dynamical fields down to one, the longitudinal wall is always 
stable. That is however an extreme choice. For typical choices of parameters, at 
cl/ct > 2 the longitudinal wall becomes unstable, and the critical bubble develops two 
kinks. Such behavior can appear with an even smaller sound speed ratio, cl/ct > V%, 
if we tune the potential to the other extreme limit. This range of sound speed ratio is 
not hard to find in real materials. 

We pick two representative forms of a (9), given by Eq. (2.12) and (2.13), to calcu- 
late the exact shapes of critical bubbles. This allows us to observe the scaling property 
of tunneling rates. When the critical bubble is deformed but still smooth, we can modify 
the standard tunneling rate formula, Eq. (1.2), through an effective Fermi velocity, 



This is quite intuitive since one particular orientation is longitudinal, and all others 
are transverse. They care about the sound speeds in their own orientations. When 
the bubble starts to develop kinks and we further increase the sound speed ratio, the 
scaling behavior changes to 



This is because the kink- development removes a large portion of the longitudinally 
oriented wall, so ci becomes less important. 

4 Analysis of the bubble shape for nucleation is identical to those of equilibrium bubbles, known as 
the Wulff construction. Earlier works [4-11] have qualitatively similar results. We further specify that 
zigzag segments of 9 C is the configuration to which the instability settles. The recognition and inter- 
pretation of the kinky bubble shape is also more transparent in our analysis. The tension previously 
studied is often expanded as a(9) = <tq + a cos 9 + b cos 29, and most analysis focused on the effect of 
a =/= 0. In our model there is a reflection symmetry — the domain-wall tension does not change when 
you look at it from the other side. Thus we always have a = 0. This makes our situation closer to a 
2D lattice model [15], where similar bubble shape was observed. 
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On top of modifying the tunneling rate estimation, our result has a practical impact. 
Typically, the experimental measurement of domain-wall tension involves measuring 
the bubble radius [16]. That is done by observing a domain-wall popping through a 
partition with holes. When it does, the radius of the hole is identified with the bubble 
radius. Our result shows that for vector fields, the orientation of that partition is 
important. Only for a longitudinally oriented partition, the popping radius can be 
identified with x given by Eq. (3.8). For other orientations, the hole and the bubble 
do not have common symmetries. Therefore the exact relation between the popping 
radius and the critical radius requires further analysis. 

We have only taken a small step toward a rich phenomenology. Given the new 
insight here, many nontrivial questions arise. How does the domain- wall move/bubble 
expand given this orientation dependence? Especially when there is a kink, can we 
expect the tip to travel at cl, leaving behind a Cherenkov-like tail of domain- walls 
bounded by cy? How does the spontaneously broken planar symmetry interact with 
impurities or other external effects? All these await future study, and may lead towards 
a more practical understanding about some exotic theories of phase transitions which 
relies on the properties of domain-walls [13,17,18]. 

Acknowledgments 

We thank Igor Aleiner, Solomon Endlich, and Erick Weinberg for stimulating discus- 
sions. This work is supported in part by the US Department of Energy, grant number 
DE-FG02-11ER41743. 

References 

[1] S. R. Coleman, V. Glaser and A. Martin, "Action Minima Among Solutions to a Class 
of Euclidean Scalar Field Equations," Commun. Math. Phys. 58, 211 (1978). 

[2] S. Coleman, "Quantum Tunneling and Negative Eigenvalues," Nucl. Phys. B298 
(1988) 178. 

[3] J. Rudnick and R. Bruinsma, "Shape of Domains in Two-Dimensional Systems: 

Virtual Singularities and a Generalized Wulff Construction," Phys. Rev. Lett. 74 (Mar, 
1995) 2491-2494. 

[4] P. Galatola and J. B. Fournier, "New Features of Two-Dimensional Soft Matter 
Domains: Dips and Quasicusps," Phys. Rev. Lett. 75 (Oct, 1995) 3297-3300. 

[5] A. J. Leggett, "A theoretical description of the new phases of liquid He-3," 
Rev.Mod.Phys. 47 (1975) 331-414. 



- 15 - 



[6] J. C. Wheatley, "Experimental properties of superfluid He-3," Rev.Mod.Phys. 47 
(1975) 415-470. 

[7] D. Vollhardt and P. Wolfle, The Superfluid Phases of Helium 3. Taylor and Francis, 
London, New York, Philadelphia, 1984. 

[8] J. B. Fournier, "Generalized Gibbs-Thomson Equation and Surface Stiffness for 
Materials with an Orientational Order Parameter," Phys. Rev. Lett. 75 (Jul, 1995) 
854-857. 

[9] J. E. Maclennan, Q. Jiang, and N. A. Clark, "Computer simulation of domain growth 
in ferroelectric liquid crystals," Phys. Rev. E 52 (Oct, 1995) 3904-3914. 

[10] J. Rudnick and K.-K. Loh, "Theory of monolayers with boundaries: Exact results and 
perturbative analysis," Phys. Rev. £60 (Sep, 1999) 3045-3062. 

[11] N. M. Silvestre, P. Patricio, and M. M. Telo da Gama, "Elliptical soft colloids in 
smectic-C films," Phys. Rev. E 74 (Aug, 2006) 021706. 

[12] A. Aguirre, M. C. Johnson, and M. Larfors, "Runaway dilatonic domain walls," Phys. 
Rev. D81 (2010) 043527, arXiv: 0911 .4342 [hep-th] . 

[13] J. T. Giblin, Jr, L. Hui, E. A. Lim, and I.-S. Yang, "How to Run Through Walls: 
Dynamics of Bubble and Soliton Collisions," Phys. Rev. D82 (2010) 045019, 
arXiv: 1005.3493 [hep-th]. 

[14] P. Ahlqvist et al, "Conifolds and Tunneling in the String Landscape," JEEP 03 (2011) 
119, arXiv: 1011.6588 [hep-th]. 

[15] I. Aleiner and R. Suris, "The Shape and Activation Energy of Critical Two-dimensional 
Nuclei on the (001) Surface of a III-V Crystal during Epitaxial Growth," Sov. Tech. 
Phys. Lett 16 (1990) 547. 

[16] D. Osheroff and M. Cross, "Interfacial Surface Energy between the Superfluid Phases 
of He-3," Phys. Rev. Lett. 38 (1977) 905-909. 

[17] R. Easther, J. Giblin, John T., L. Hui, and E. A. Lim, "A New Mechanism for Bubble 
Nucleation: Classical Transitions," arXiv : 0907 . 3234 [hep-th] . 

[18] I.-S. Yang, S.-H. Tye, and B. Shlaer, "Classical Transitions in Superfluid Helium 
3," ArXiv e-prints (Oct., 2011) , arXiv: 1110.2045 [cond-mat . other] . 



- 16 - 



A. Examples of different potentials. 



Here we present some analytical approaches and the numerical evaluation of the domain- 
wall tension given by the relaxation method [12-14]. The first potential we study is a 
double well in one direction and a quadratic in the other direction. 

V{K <Py) = Aftfy + + \P€ ■ (A.l) 

This is qualitatively similar to Eq. (2.2) when m 2 > 0. The differences are some 4th 
order terms involving (j) x , which is not very important when (3 > stabilizes a trajectory 
near <p x = 0. The two potentials can be roughly related by 

,2 ™2 

+ b\H\ 



2 2 

^ - ^ • (A.2) 
2 2 v ; 



The two degenerate minima sit at (0, ±y ^-). For the purely longitudinal (or trans- 
verse) wall oriented along the y (or x), we can solve the problem analytically and get 
the exact value of the tension. 
Longitudinal wall (8 = 0): 



^(*,y) = ^taiihl-^J (A.3) 

Mx,y) = o (A.4) 

2V2fu 3 c L 

° = a ^ = — 3A — ' ^ ^ 



Transverse wall 



^(z,,) = ^tanh(-^J (A.6) 
<j) x (x,y)=0 (A.7) 

a = g{<k/2) = v ^ • (A.8) 

For other orientations of the wall, we can evaluate the tension numerically. Before 
that, we can analyze two extreme cases. Using the method of rotating the potential as 
described in Sec. 2, we have 

r 2 r 2 

U T i/2 , U L i/2 



a(9) = I dx + -±<ft + V 9 (4> x , <j> y ) 



f c 2 c 2 

/ dx + y <fiy + v {4>x cos 9 + (j) y sin 9, <p y cos 9 - <j> x sin 9) . (A. 9) 
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When j3 — > oo, we effectively have a single field problem with 

*=^ = ^> (A.10) 

cos 9 sin v 

such that 

f , c r cos 2 # + sin 2 9 - n rr/ - . . . 

<r(0) = / rfa; ^ ^ 0' 2 + V(0, 0) . (A.ll) 

Clearly, this gives us Eq. (2.13). 

The other extreme limit is /3 — > 0, at which the potential is flat in the X direction. 
The two degenerate vacua approach two separated lines. Moving along these lines 
contributes nothing to the tension. As shown in Fig. 5, the path that minimizes the 
tension involves first moving along these lines to an appropriate angle 0, then connecting 
through a straight line. The tension of this path is a function of both 9 and through 
the orientation dependence in Eq. (2.13), and a simple projection of the length. 

<r(0) 1 



a(9, 0) = Cl cos , g _f \/<?L cos 2 + 4 sin 2 . (A.12) 



Minimizing this with 0, we have 



c|. cos 9 



<t>m{9) = arccos . (A. 13) 

y Cy cos 2 9 + c\ sin 9 

So the tension in this case should be 

o{9) = a\9^ m {9)\ . (A.14) 

We can apply the analysis in Sec. 2. 2 and calculate the stability condition for the 
flat longitudinal wall. 

We can see that the wall becomes unstable as soon as cl > V^ct- 

We next provide several plots with the numerical values on top of the three possible 
fits, Eq. (2.12), (2.13) and (A.14). Fig. 6 shows that the two extreme limits indeed 
fit very well with our analysis. Fig .(7) shows that with a more moderate choice of 
parameters, Eq. (2.12) is quite reliable independent in various sound speeds. 

In the end, we provide a much more complicated potential as in Fig .(8). It has 
a general slope in the X direction and two minima located at (0.001, ±2.498). So 
trivially, the interpolation path will always involve both fields. 



V(<p x ,<p y ) = e^ \l-Sexp 



,05,-7*1. 

sinf— 

7*2 — 7*1 ' 



tanh 2 ^v-^v-^) 

(A.16) 
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Figure 5: The cl and ct axes are the directions in which the field has purely longitudinal 
and transverse sound speeds. The two dots represent the two discrete vacua. In the limit 
/3 — > 0, the dashed lines through them are almost in the vacuum, too. The important portion 
of the domain-wall is the red (thick) path from one line to the other, which is free to pick the 
best orientation <fi. 



c L = 15.8 c T =5 A = 1 n =4 /3 = 0.01 c L = 15.8 c T =5 A = 1 fj -4 fj = 500 

cr(6) 




Figure 6: The numerically calculated values of the tension for a double- well potential are 
shown in dots. The three analytical fits: Eq .(2.12) is the dashed line, Eq. (2.13) is the 
dot-dashed (blue) line, and Eq. (A. 14) is the solid (red) line. We can see that in the for /3, 
Eq. (2.13) is a good fit, and for small /3, Eq. (A. 14) is a good fit. 

We numerically evaluated the tension for various orientations and plotted it against 
the three analytical fits in Fig. (9). The overall shape can be quite different from any 
equation given in this paper. In particular, note that in the right portion of Fig. (9) 
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c L = 4V 2 c r =4 A = 1 /i=4 /? = 5 



c L = 10.2 c T =2 A = 1 yu =4 j8 = 5 



CT(0) 

341 




Figure 7: The numerically calculated values of the tension for a double-well potential are 
shown in dots. Again the three analytical fits: Eq .(2.12) is the dashed line, Eq. (2.13) is the 
dot-dashed (blue) line, and Eq. (A. 14) is the solid (red) line. The two figures use the same 
potential but different sound speed ratios. 




Figure 8: The more complicated potential introduced in Eq .(A. 16), with q = 0.5 , n = 
—2.5 , T2 = 2.5 . We will use two values of S, 1.1 and 0.9, but that makes no visual difference. 

the longitudinal domain- wall (actually an open set near 8 = 0) does not exist. 5 



5 This comes from the same reason as described in [12]. The interpolation path breaks into two 
parts, connecting each vacuum individually with the — <j) x region. For vector fields, such runaway 
behavior also acquires an orientation dependence. 
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Figure 9: From the potential given by Eq .(A. 16), we again compare the numerical <j{6) 
with the three equations. In the left figure we have 5 = 1.1. In the right figure we have 5 = 1 
and for some orientations the domain-wall does not exist because the path runs away toward 
the — 4> x direction. 
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